Neighbor-junction state effect on the fluxon motion in a Josephson stack 
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We study experimentally and theoretically the influence of phase-whirling (resistive) state in one 
junction of a two-fold Josephson stack on the fluxon motion in the other junction. In experiment, 
we measure the fluxon velocity versus current in one junction as a function of the state (Meissner 
or resistive) of the neighboring junction. The analysis, made for the limit of high fluxon density, 
shows that the interaction with the resistive state results in an increase of the effective damping 
for the moving fluxon and, therefore, in reduction of its velocity. Numerical simulations confirm 
this result for various fluxon densities. The experimental data are in good agreement with the 
theoretical predictions. In addition, the fluxon step measured experimentally has a rather peculiar 
structure with back and forth bending regions which is understood as a manifestation of the photon 
absorption in the neighboring junction. 
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I. INTRODUCTION 

Stacked long Josephson junctions (LJJ's) have re- 
cently received much attention iSjiwse they show a vari- 
ety of new physical phenomenaErO in comparison with 
single LJJ's and have potential for applications as a nar- 
row linewidth powerful oscillators for mm and sub-mm 
wavebandsQ. The naturally layered high-T^ superconduc- 
tors (HTS) cajp. be described as intrinsic stacks of Joseph- 
son junctionsQ. Therefore, study of fluxon dynamics in 
artificial stacks can help to understand the phenomena 
which take place in HTS. 

The inductive coupling model describing the dynam- 
ics of Josephson phases inpJV inductively coupled LJJ's 
was derived by Sakai et alu. Experimental investigation 
of stacked junctions became possible aftra; the progress 
achieved in (Nb-Al-AlOx)Ar-Nb technologyB which, at the 
present stage, allows to fabricate stacks with up to about 
30 Josephson tunnel junctions having parameter spread 
between them of less than 10%lI. Initially, the interest 
was concentrated on investigation of the simplest sym- 
metric fluxon states since they are promising for oscilla- 
tor applications. Later on, it was found that it is very 
interesting to understand the asymmetric stateSntjecause 
they show rather nontrivial nonlinear dynamicgj'El. Such 
asymmetric states are also of practical importance, be- 
cause multilayered oscillators most probably will operate 
in a regime when only the majority but not all of the 
junctions oscillate coherently while the other juiictions 
are in resistive or opt synchronized flux-flow statoj. 

In a recent workQ it has been shown that the dynamic 
state of one junction in a 2-fold stack affects the static 
properties of the other junction. As a next step, it is 
interesting and important to understand how different 



dynamic states in one LJ J affect the dynamics of fluxons 
in the other LJJ. In particular, the goal of this work is 
to study the dynamics of a fluxon in one LJJ when the 
neighboring LJJ is in the resistive state and compare it 
with the case when the neighboring LJJ is in Meissner 
state. We call the resistive state a "phase-whirling" state 
because the Josephson phase difference rotates very fast 
and nearly uniformly, in rough approximation. Such a 
dynamic state often occurs in experiments and, there- 
fore, it is important to understand and describe it ad- 
equately. In fact, in the early experiments with stacks 
it was somewhat naively supposed that the voltage of 
Flux-Flow Step (FFS) in one LJJ does not depend on 
the state (Meissner or resistive) of the other LJJ. In fact 
this is true only for the asymptotic voltage of FFS. Here 
we show that in the presence of the "phase-whirling" so- 
lution in one of the junctions the actual flux-flow voltage 
across the other LJJ gets lower. 

In the section II we present the experimental data 
which clearly show that in a 2-fold stack the switching 
of one junction from the Meissner state to the phase- 
whirling state decreases the velocity of a fluxon moving 
in the other junction and, therefore, the dc voltage across 
it. The analytical approach which explains the observed 
decrease of fluxon velocity in the limit of high fluxon den- 
sity is developed in section III. The results of numerical 
simulations confirm analytical results and are shown in 
section IV. The results of the work are summarized in 
section V. 



II. EXPERIMENT 

In order to investigate the influence of the phase- 
whirhng state in the neighboring junction on the fluxon 
dynamics, we have chosen the most clean ring-shape (an- 
nular) LJJ stack geometry. Due to magnetic flux quanti- 
zation in a superconducting ring, the number of fluxons 
initially trapped in each annular junction of the stack is 
conserved. The fluxon dynamics can be studied here un- 
der periodic boundary conditions which exclude possible 
complicated interference of the fluxon with the junction 
edges. 

Experiments have been performed with 3 different (Nb- 
Al-A102;)2-Nb annular LJJ stacks prepared in 2 different 
technological runs (2 samples in one run and the third 
sample in another run). The sample geometry is shown 
in Fig. |l|. Two annular LJJ's are stacked one on top of 
the other, with bias leads attached to the top and bottom 
electrodes. The physical parameters of all samples, mea- 
sured at T = 4.2 K, are summarized in Tab. |. The stacks 
were designed with extra contacts to the middle super- 
conducting electrodeEEl so that the voltages across each 
LJJ can be measured separately. The inner diameter of 
all stacks was D — 122.5 /im and the width W — 10 /im. 
Due to technological difhcultics of making a stack of iden- 
tical LJJ's with contacts to the middle electrode, the 
two stacked junctions had rather substantial difference 
in quasiparticle (subgap) resistance -Rqp. The normal- 
ized circumference of the ring was t:D / \j = L/Xj « 15, 
where \j is the Josephson penetration depth, which was 
approximately equal in both junctions. Measurements 
were performed in the temperature range 4.2-5.8 K. 

In stacked annular LJJ's, clean trapping of a single 
fluxon in a desired junction is rather difficult due to the 
asymmetry of the required state [1|0]. In the particular 
case of the 3 samples mentioned above, the asymmetry 
in the junction's resistance allowed to trap the fluxon in 
the desired [1|0] state without many efforts just by ap- 
plying a small bias current through one of the junctions 
during cooling the sample below the critical temperature 
Tc- After every trapping attempt, the resulting state was 
checked. The I-V characteristic (IVC) of both LJJ's 
were traced simultaneously in such a way that the cur- 
rent was applied through the whole structure (through 
two junctions connected in series) and the voltages were 
measured individually across each LJJ. The wanted state 
[1|0] with a fluxon in one junction and no fluxon in the 
other junction was identified by a simultaneous observa- 
tion of a small critical current Ic and fluxon step with 
the smallest asymptotic voltage '^ 20 /xV in LJJ"^, and a 
large critical current in LJJ^. Both the current ampli- 
tude of the fluxon step I^^^{H) and the critical current 
J^ {H) are expected to have their maxima at zero applied 
magnetic field H — 0. To check that we have clean fluxon 
trapping, i.e. that the fluxon is trapped in a LJJ and not 
accompanied by the parasitic Abrikosov vortices in the 
superconducting films surrounding LJJ, we checked the 



dependences I^'^{H) and Imax{H) after each trapping 
attempt and repeated it until these dependences were 
symmetric. 

The main experimental result of the paper is shown in 
Fig. 0. It is IVC's of both LJJ's of the sample #2 traced 
at T « 5 K using rather complex current sweep sequence. 
Note, that the voltage scales of two IVC's in Fig. are 
different and shown on the bottom axis for V and on 
the top axis for V^ . The sweep starts at the bias point 
A where J = and V = and a fluxon is trapped in 
LJJ^ (state [1|0]). When the current is increased up to 
about / — 0.69 mA (point B in Fig. g), LJJ^ switches to 
the fluxon step, while LJJ^ still remains in the Meissner 
state. Ideally, the LJJ"^ should switch to the fluxon step 
at zero bias current since any non zero current applied 
should drive the fluxon around the stack. In our case the 
fluxon is pinned, most probably near one of the contacts 
to the middle electrode, and only current / — 0.69 mA 
can tear it away from the pinning center. With the fur- 
ther increase of the bias current the LJJ^ follows the 
fluxon step which corresponds to the fluxon rotating in 
the ring, and the voltage across LJJ"* is proportional to 
the fluxon rotation frequency according to the Josephson 
relation. 

In an ideal single annular LJJ, the fluxon step has a 
relativistic nature and its slope approaches infinity when 
the fluxon moves with velocity u close to the Swihart ve- 
locity cq. In the stack with different jc or with inhomo- 
geneities, the fluxon's velocity can exceed the Swihart ve- 
locity. This results in the emission of the electromagnetic 
waves travelling behind moving fluxon (Cherenkov rjadi- 
ation) and in a finite slope of the step at any velocityuEiJ. 
If the length of the emitted radiation tail is comparable 
with the circumference of the LJJ, the resonant struc- 
tures (small steps) can appear on the top of the fluxon 
step. Such steps are visible on the top of the fluxon step 
in Fig. g and are outlined by the circle. At / « 2.67 mA 
(point C in Fig. 0) both junctions simultaneously switch 
to the resistive state (gap voltage) with rapidly whirling 
Josephson phase. Such a simultaneous switching is called 
current locking. We studied it in detail for stacks of linear 
geometry in Ref. @. Analyzing the dependence of critical 
current I^ (H) and maximum current of the fluxon step 
-^inax(^) on magnetic field, we conclude that the current 
locking was driven (initiated) by LJJ"* (if LJJ^ is kept 
in the resistive state, I^ is substantially higher). 

When both LJJ's are in the resistive state, we inverse 
the direction of the sweep, i.e. start reducing the bias 
current. At / = 1.08mA (point D in Fig. ||), LJJ"* 
switches from the resistive state to the fluxon step while 
LJJ^ still stays in the phase-whirling state. We denote 
such state of the stack as [l|i?]. The voltage V^ in the 
[1\R] is by about 16% smaller than V^ in the [1|0] state 
at the same bias. In fact, this is one of key observations 
in our study. 

At this point there are two possibilities: first, continue 
to decrease the bias current down to zero or, second, 
increase the current and trace up the single fluxon step 



for [1\R] state. 

If we continue decreasing the bias current, at / = 
0.916 mA (point E in Fig. ||) LJJ^ switches from the 
resistive state (McCumber branch) to the Meissner state 
and the overall state of the stack becomes [1|0]. This 
causes the voltage V^ to increase and become equal to 
voltage of the fluxon step which we traced in the begin- 
ning of the bias current sweep. The fact that switching 
of LJJ^ caused a voltage jump across LJJ^ is marked 
in Fig. H by dotted arrow. Thus, we demonstrated ex- 
perimentally that the change in the state of LJJ^ af- 
fects the velocity of fluxon moving in LJJ^. Further de- 
crease of the bias current results in the fluxon pinning at 
/ — 0.470 mA (point F in Fig. g) and in the zero voltage 
across both LJJ's. 

The second possibility is, being in the bias point D, to 
increase bias current and trace the fluxon step of the [1\R\ 
state up. This step has rather peculiar shape as shown in 
Fig. g. In addition to the common trend to have smaller 
voltage than the fluxon step in the [1|0], the step in [l\R] 
state bends back and forth which possibly implies some 
interesting physics behind it. We also noticed that as we 
trace both IVC's in [R\R] state from bias point C down 
to bias point D and then in [1\R] state from point D 
up to bias point G, the voltage across the LJJ^ has a 
small hysteresis at voltages equal to the sum of the gap 
voltages of the superconducting electrodes constituiting 
the LJJ^. This small hysteresis is shown magnified in 
the inset of Fig. 0. The voltage across LJ J-^ is somewhat 
smaller in the [1\R] state than in the [R\R] state. As soon 
as LJJ'^ switches to the resistive state (point G in Fig. |^ 
and dotted arrow in the inset), this difference vanishes. 

We propose the following explanation for the observed 
back bending. As we increase current starting from point 
D up to / w 1.5 mA, corresponding to nearly vertical 
slope of the fluxon step, the voltage V^ increases from 
55 /iV up to 1.9 mV, i.e. approaches the gap voltage. 
The fluxon motion in LJJ , due to the coupling between 
the junctions, causes oscillations of Josephson phase and, 
therefore, of electric and magnetic fields in LJJ^ JChis 
leads to photon-assisted tunneling (PAT) effectEj£j in 
LJJ-^. This effect was earlier observed by GiaeveitJ also 
using a stack of two junctions. The characteristic fre- 
quency oj of photons absorbed in LJJ^ is equal to the 
fluxon rotation frequency in LJJ . Due to PAT, one ex- 
pects to observe a step at the gap sum voltage decreased 
by huj/e = 2V^. In the low bias region shown in the 
inset of Fig. |2|the gap sum step is not well defined which 
results in somewhat weaker gap suppression result. In 
fact, the maximum suppression of the voltage V^ we 
have found is about 20 /iV that corresponds to the top 
of the back bending region at / w 2.6 mA. Since the 
Josephson voltage in LJ J^ is also about 20 /iV, its effect 
on V^ is only about 50% of the expected PAT step volt- 
age change. Thus, as we increase current, the gap voltage 
in LJJ^ decreases due to photon-assisted tunneling. The 
resulting decrease of fluxon step voltage of LJ J^ is asso- 
ciated with the appearance of an additional dissipation 



channel due to PAT. Since the PAT step on IVC is lim- 
ited in voltage (by 2V ) and in current amplitude (oc to 
the amplitude of the first Josephson harmonic which, in 
the case of fluxon motion, saturates at some bias), the 
bending to the right caused by Cherenkov radiation ap- 
pears to be stronger at / > 2.6 mA so that the fluxon 
step of LJJ^ gains the positive slope again. The differ- 
ential resistance at the top the fluxon step in [l\R] state 
is rather high and no resonances are observed. This pic- 
ture is typical for fluxon with Cherenkov radiation tail 
moving in a media with high dissipation. 

The negative bias part of the IVC reproduces all the 
features described above for the positive half except for 
the small hysteresis between bias point B and F. Very 
similar IVC's were found for other 2 measured samples. 
The minor difference was in the particular values of the 
bias current in the points B, C, D, E, F and G which 
were also dependent on T. All samples showed that the 
voltage V^ of the fluxon step in the [1|0] state is some- 
what higher than the voltage of the same step at the 
same bias in the [1\R] state. The fluxon step in the [1\R] 
state showed back and forth bending for all samples. The 
hysteresis between points B and F was intersecting with 
the hysteresis between points D and E for some samples 
and temperatures, so we had to use even more complex 
sweep sequence in order to trace out all possible dynam- 
ical states. 



III. THEORY 

The main objective of this section is to analyze the ori- 
gin of the decrease in the fluxon velocity in one LJJ due 
to the switching of the neighboring junction into the re- 
sistive state. Here we use the standard RSJ model which 
does not take into account the dependence of the dissi- 
pation on voltage a.t V ^ Vg ^ 2.4 mV. Thus the gap 
related effects like PAT discussed above are neglected. 

The fluxon dynamics in the system under investigation 
can be descdbed in the framework of the inductive cou- 
pling modelEI which for the case of two coupled junctions 
takes the form: 
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where </> and -0 are the Josephson phases across the 
LJJ"^'-^, respectively, — l<S'<Oisa dimensionless 
coupling parameter, a <C 1 is the damping coefficient de- 
scribing the dissipation in the system due to quasiparticle 
tunneling, and 7 — j / jc is the normalized the density of 
the bias current flowing through the stack. The coordi- 
nate X and time t are measured, respectively, in units of 
the Josephson length \j and inverse plasma frequency 
LOp"^ of single-layer LJJ. Most of relevant parameters of 
the junctions, such as effective magnetic thicknesses and 



specific capacitances, for the sake of simplicity are as- 
sumed to be equal in both LJJ's. 

To understand the origin of an additional friction force 
we start from unperturbed (without r.h.s.) Eqs. (|l|) and 
(0) and use the force balance equations to derive the 
shape of the IVC. We do not directly solve Eqs. (|l|) and 
(0), but, rather, use trial solutions for 4>{x, t) and ipi^i t)i 
which approximate the Josephson phase profiles in the 
states [1|0] or [1\R]. The choice of the trial functions 
is suggested by the results of numerical simulations pre- 
sented in the following section. To simplify the mathe- 
matics and concentrate attention on the physical sense, 
the dense fluxon chain approximation is used. For this 
case we adopt the following trial solutions: 

0(x, t) = H{x - ut) + Ar.m siu [if (a; - ut)] ; (3) 
^''\x,t)^B^s\n[H{x~ut)]- (4) 

V'^'fa;, t) ^ Lot ^ sin(wi) -I- Br sin [Hix — ut)] , (5) 

where -0™ and ip^ are phases in LJJ^ in the Meissner 
state and resistive state, accordingly (these are the two 
cases which we are going to compare); u is the veloc- 
ity of the fluxon chain; Ar^rm Br^m *^ 1 are the con- 
stants which we are going to determine for resistive and 
Meissner states, respectively; H is the average normal- 
ized magnetic field in the LJJ^: 



H = 



2nN 



(6) 



where N is the number of fluxons trapped in annular 
LJJ-^ (for [1|0] and [1\R] states TV = 1), and £ = L/Xj 
is the normalized length of the junctions. Dense fluxon 
chain approximation implies that H ^ 1. 

Substituting Eqs. (|)-(|) into Eqs. @) and @, and 
using the following approximations (which are justified 
in our case) 



sin (j) « sin [H{x — ut)] ; 

sinip"" Ki B sin [H{x-ut)]; 

sinip^ « sin(ajt). 



(7) 
(8) 
(9) 



we arrive to the equations from which we can determine 
Ar^m and Br.m- The final result for the [N\R] (resistive) 
state is 



Br — 



D 

-S 



where we introduced notations 

D = \-u^{l-S^)] 
Q={l-u^)''-u^S^. 

Note, that both D > Q and Q > for u < c^ 
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The result for the [A^|0] (Meissner) state is 
DH'^ + 1-52 

Am — ^ 
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(14) 
(15) 



To calculate IVC, 7(u), we write the force balance 
equation 



2nN^ = F^ + Ff . 



(16) 



Here F^'^ are the friction forces which develop in 
LJJ'^'-^. The expression for the friction force is well 
known from the perturbation theoryEj 



Fa = a 



dx, 



(17) 



where we will use i = 2ttN/H following from (||). Since 
we are interested in the average friction force to get an 
IVC, we have to perform friction force averaging in time: 



T 



Fa{t)dt. 



(18) 



Since there are two characteristic frequencies in the 
system, fluxon (Josephson) frequency and the phase- 
whirling frequency of the resistive state, we have to 
choose the averaging interval T in Eq. (08) so that it 
will contain an integer number of periods of each fre- 
quency i.e. T = 2-Kklu} = 2-Km/Hu, and fc, m being the 
integer constants. After averaging, we get the following 
expressions for friction forces 



F^ = TiNaHu {Al„, + 2) ; 
F^ ^TiNaHuBl^. 



(19) 
(20) 



All information about the actual state is contained in 
Jr m calculated above for the Meissner and 



J\.r m ancl iJr 



resistive state of the LJJ'^. 

Finally, we insert (n9p and ( |20| ) into the force balance 
equation (nq) and get IVC's 



7r,m(") 



aHu 



K, 



Bi 



(21) 



Now we can prove that IVC for the [A^ji?] state is shifted 
to the region of lower velocities in comparison with the 
IVC for [A^|0] state, i.e. that 



5{u) — ^r{u) — ^m(u) > for all \u\ < c_. 



(22) 



Substituting Ar,m and Brm from Eqs. (|10|), (yj), (|lj), 
(|l5| ) into the expression (£l|) and using the obtained ex- 
pressions for jr,m{u) in (|22|) we get 



S{u) 



au (XiH'^ + X2) 

2ipQ^~{ipQnyf 



(23) 



where Xi and X2 are defined as 

Xi ^ 2QD [S^ -Q{l-S^)+ L»2] 



(24) 



X2 = D^(D^ + S^)-(l-S^)^Q^, (25) 



Obviously, Eq. (|2^) is positive when both Xi and X2 
arc positive. To prove the latter, we express D^ as a 
function of Q using Eqs. (^2|) and ( [l3| ) 



D^ = Q (1 _ 5-2 



5^ 



Substituting (||) into (|2|) and (|5|) we get 

Xi = AQDS^ > 0; 

X2 = S^ [3Q (1 - 52) + 2S'2] > 0. 



(26) 



(27) 
(28) 



Thus ( p^ ) is proved and 7r('«) > 7m (u) for any u < c_. 

This result is in agreement with our experiment and 
simulation (see the following section). From physical 
point of view, the origin of the effect lays in the difference 
between Eqs. (g) and (|9|) where the main term depends 
on the state of LJJ'^ resulting in different phase profile 
and different friction force for [1|0] and [1\R] states. The 
IVC's of fluxon steps ^"^(7) for the states [llO] and [1\R] 
calculated using Eq. (|2^) are shown in Fig. ^. According 
to the calculations presented above the difference 6{u) 
diverges as u approaches c_. Actually, for this case of 
a single fluxon in our relatively long junction, the dense 
fluxon chain approximation is not fully valid so one needs 
to perform more exact analysis or numerical simulations. 
The region of validity of our approximation is A^ < 1 
and Br < I. Using Eqs. (|l(]| ) and ( |ll] ) for the same pa- 
rameters as in Fig. we get that our approximation is 
valid up to w w 0.809 while c_ « 0.816. 



IV. SIMULATION 

To check the limitations of the analysis presented 
above we performed direct numerical simulations. Our 
simulations show that the effect observed in experiment 
and explained in the framework of high fluxon density 
approximation exists for any fluxon densities and any 
velocity even very close to c_. Another advantage of 
the simulation over the analytical approach is that the 
simulation fully reproduce the dynamics of the inductive 
coupling model and, therefore, all the effects possible in 
its framework. 

The numerical procedure works as follows. For a given 
set of LJJ's parameters we simulate the IVC of the sys- 
tem, i.e. calculate V^{j) and ^^^^(7) while increasing 7 
from zero up to 1. To calculate the voltages V^{'-f) and 
^^(7) for each value of 7, we simulate the dynamics of 
the phases (j)"^'^{x, t) by solving the Eqs. (fl^) and (0) with 
the periodic boundary conditions: 



0^^^(O,i)=0^^^(£,t) + 27rAr 



tA,B 



i.A,B 



A,B I 



C'"(O,i)=0,'"(^,t) 



(29) 
(30) 



numerically using an explicit method [expressing 
(|)A,B(^^ + Ai) as a function of (j)^^^{t) and (j)^^^{t - Ai)], 
and treating (l)xx with a five point, <^tt and (^t with a three 
point symmetric finite difference scheme. Numerical sta- 
bility was checked by doubling the spatial and temporal 
discretization steps Ax and Ai and checking its infiuence 
on the fiuxon profiles and on the IVC. The discretization 
values used for simulation were Ax = 0.01, Ai = 0.0025. 
After simulation of the phase dynamics for Tq = 20 time 
units we calculate the average dc voltages V^^^ during 
this time interval as 



yA,B ^ 1 
T 



,A.B 



(/)f'^(i)di = 



(/)^'^(r)-0^'^(o) 

T 



(31) 



For faster convergence, we use the fact that V"^'^ does 
not depend on x and, therefore, we also take advantage 
of the spacial averaging of the phases (/)^'^ i n (pl| ). 

When the values of V^^^ are found from (|3l|), the dy- 
namics of the phases (f> '^{x,t) is simulated further dur- 
ing 1.2ro time units, the dc voltages V"^'^ are calculated 
for this new time interval and are compared with the 
previously calculated values. We repeat such iterations 
further increasing the time interval by a factor 1.2 until 
the difference in dc voltages |y(1.2"+i T) - ^(1.2" T)\ 
obtained in two subsequent iterations becomes less than 
a given accuracy SV = 10"'^. The particular value of 
the factor 1.2 was found to be quite optimal to provide 
fast convergence as well as more effective averaging of 
low harmonics on subsequent steps. Very small value of 
this factor, e.g. 1.01 can result in very slow convergence 
in the case when (/)(i) contains harmonics with the pe- 
riod comparable or larger than T. Large values of the 
factor, e.g. 2 or higher, will consume a lot of CPU time 
already during the second or third iteration even when 
the convergence is good. After the voltage averaging for 
current 7 is complete, 7 is changed by a small amount 
S^ to calculate the voltages in the next point of the IVC. 
As initial conditions here we use a distribution of phases 
(and their derivatives) achieved in the previous point of 
the IVC. 

An example of calculated IVC is shown in Fig. H To 
trace both the Meissner and resistive states we use the 
following sweep sequence: 7 increases from up to 1 with 
a step ^7 = 0.01, then decreases down to 0.5 with a step 
(57 — 0.01, and further down with a step 5'^ = 0.002 until 
the state [-/V|i?] is reached as shown in Fig. Q. From this 
point, we either sweep further down to 7 = or up to 
7=1 until both junctions switch to the resistive state. 
The decrease of voltage in [1\R] state in comparison with 
[1|0] state is very clearly seen in Fig. H. We performed 
this kind of simulations for the wide range of the param- 
eters i.e. for IS"! = 0.1, 0.2... 0.5, and TV = 1, 2 . . . 5 
(in total 25 pairs of IVC's) and found similar IVC's in 
all cases. For relatively dense fluxon chain N/L = 1 we 
also compared the IVC's obtained by means of numeri- 
cal simulation with IVC's derived analytically and found 
a good agreement as shown in Fig. 0. Small difference 



in the slope of analytical and numerical IVC's is related 
to final density of fluxons in simulation H = 2t: while 
theoretical curve corresponds io H ^ 1. 

The limitations of our model due to voltage indepen- 
dent loss term a prevented proper calculation of the up- 
per part of the fluxon step in [1\R\ state. Here numerical 
curve shows a series of small voltage jumps (see Fig. ^) in 
LJJ that are not observed in experiment. These jumps 
are related to excitation of fluxon-antifluxon states in 
LJJ-^ at high voltages. Such states were not found in ex- 
periment, most probably due to the increased dumping 
in LJJ^ at gap voltage. 



V. CONCLUSION 

We investigated experimentally and theoretically the 
motion of a fluxon in one of two magnetically coupled 
long Joscphson junctions. Two different cases are stud- 
ied: (1) when the neighboring junction (one which does 
not contain any fluxon) is in the Meissner state, and (2) 
when the neighboring junction is in a phase-whirling (re- 
sistive) state. We found that the phase-whirling state 
in LJJ^ slows down the fluxon motion in LJJ"^ and re- 
sults in a shift of the fluxon step to lower voltage. This 
effects is detected experimentally, reproduced in simula- 
tions based on the inductive coupling model, and derived 
analytically in the high fluxon density approximation. In 
addition, experiment shows quite peculiar back and forth 
bending of the fluxon step in [l|i?] state which we explain 
as a result of increased dumping due to photon assisted 
tunneling effect in LJJ^. The results of our study are 
also relevant for characterization of stacked Josephson 
junctions with large number of layers. 
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FIG. 1. Two coupled (stacked) long Josephson junctions 
of annular geometry with one fluxon trapped in the top junc- 
tion. For symmetry, we used two voltage probes attached to 
the middle superconducting electrode. Dimensions are not to 
scale. 



FIG. 2. Experimentally measured IVC of LJJ containing 
a fluxon (open circles) and LJJ^ (solid dots). The IVC's 
are plotted with different voltage scales: top axis shows the 
voltage across LJJ^ which bottom axis — across LJJ . This 
plot also shows a magnified view of the IVC corresponding to 
LJJ^ at T/«2.2mV. 
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FIG. 3. IVC's of mo] and [l\R] states as predicted by 
analytical model Eqs. (Glh (shown by lines) and obtained as a 
result of numerical simulation (symbols) for S = —0.5, N — 5, 
£ = 5 {H = 2-k). 



FIG. 4. IVC's obtained by means of numerical simulation 
for £ = 5, a = 0.05, N ^1, S = -0.5. The arrows show the 
directions of the sweep. 
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sample 


#1 


#2 


#3 


V, (mV) 


2.38/2.54 


2.37/2.54 


2.51/2.61 


Ic (mA) 


~6.0 


-6.2 


~7 


^QP (fi) 


0.6/5.1 


0.4/4.0 


1.1/5.4 


Rn (fi) 


0.19/0.18 


0.16/0.17 


0.18/0.18 


A/g (mA) 


10/12 


9/13 


10/10 



TABLE I. Physical characteristics of stacks measures at 
T — 4.2 K. Two numbers separated by slash, are related to 
top and bottom LJJ of the stack, respectively. 
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